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Abstract -We study topological properties of phase transition points of one-dimensional topo¬ 
logical quantum phase transitions by assigning winding numbers defined on closed circles around 
the gap closing points in the parameter space of momentum and a transition driving parameter, 
which overcomes the problem of ill definition of winding numbers on the transition points. By 
applying our scheme to the extended Kitaev model and extended Su-Schrieffer-Heeger model, we 
demonstrate that the topological phase transition can be well characterized by winding numbers of 
transition points, which reflect the change of the winding number of topologically different phases 
across the phase transition points. 
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Introduction. — Due to the rapid progress in the 
study of topological insulators and superconductors mu, 
theoretical studies of topological phases and phase transi¬ 
tions have attracted great interest mm- Unlike conven¬ 
tional quantum phase transitions (QPTs) [12], topological 
QPTs do not accompany symmetry breaking, but involve 
the change of ground-state topological properties. Differ¬ 
ent topological states can be classified by quantized topo¬ 
logical invariants, while the different phases in conven¬ 
tional QPTs are distinguished by order parameters, which 
generally take continuous values. Conventional QPTs can 
be classified by singularity properties of ground-state en¬ 
ergy at the phase transition point, but this method does 
not unveil the topological properties of the phase tran¬ 
sition point, and can not distinguish a topological phase 
transition from a conventional one. On the other hand, 
topological invariants are usually defined for quantum 
states protected by nonzero energy gaps, hence they are 
ill defined at topological phase transition points, where 
the band gaps close at some specific points in the Bril- 
louin zone (BZ). Nevertheless, it was recently indicated 
that the topological properties of phase transition points 
can be characterized by introducing topological invariants 
defined on closed curves surrounding the transition points 
in the enlarged parameter space d. which provides a 


scheme to study the topological properties of phase tran¬ 
sition points. 

Since its discovery d> the Berry phase has played an 
important role in many quantum phenomena m- A well- 
known example is that the quantized transversal conduc¬ 
tivity in the quantum Hall effect can be measured by the 
Chern number, i.e., the integral of the Berry curvature 
over the two-dimensional BZ. The topological properties 
of one-dimensional (ID) systems can be classified by the 
Berry phase across the BZ (i.e. the Zak phase HZ]), which 
can be experimentally measured in ID optical lattices [18] . 
The Zak phase 7 takes the value of n with a modulus of 2n 
for topologically nontrivial systems which hold a pair of 
degenerate edge states mu, while there are no edge states 
in topologically trivial systems with 7 = 0. However, for 
Z-type ID topological systems, there may exist different 
topological phases with more than one pair of degenerate 
edge states [715] . and can not be fully classified by the Zak 
phase. For ID topological systems, the winding number is 
a more convenient quantity to characterize the topological 
properties of the Z-type systems, but it also suffers from 
the problem of ill definition on topological phase transition 
points. 

In this letter, we first give a definition of the winding 
number for the phase transition point in the enlarged pa- 
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rameter space by treating the phase driving parameter as 
an additional parameter besides the momentum. Such 
a definition overcomes the ill-definition problem at the 
transition point as it is defined on a closed detour path 
surrounding the gapless transition point and enables us 
to study the topological property of the transition point. 
Via some straightforward arguments, we find a simple rela¬ 
tion between the defined winding number for the transition 
point and the winding numbers on two sides of the tran¬ 
sition point defined in the BZ, which explains why it can 
be used to classify various topological phase transitions. 
To exemplify our theory, we then apply our scheme to 
study the extended Kitaev’s p-wave superconductor model 
and the extended Su-Schrieffer-Heeger (SSH) model which 
both exhibit rich phase diagrams. 

Winding number for phase transition point. — 

We consider a general two-band ID Z- type topological 
system, whose Hamiltonian in momentum space only con¬ 
tains two of the three Pauli’s matrices, and can be written 
in the form of 


— h()I T h x &x T dyCTy (f) 


after some rotations, with I the unit matrix and a the 
Pauli’s matrices. Here 77 represents a phase transition driv¬ 
ing parameter, k is the ID momentum, and ho, h x , h y are 
generally functions of parameters 77 and k. The winding 
number for ID systems can be associated to the Zak phase 


n2n 

7 = / dk(u k \id k \u k ), 
Jo 


( 2 ) 


where u k denote the occupied Bloch states of the Hamil¬ 
tonian. The eigenstate for the lower band then has the 
= , — 1 ) T , thus the Zak phase 


form of u k = -^= I 
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where 


v = 


1 

2 n 


t^X dhy /iy dh X 

hi + h\ 


(3) 

(4) 

(5) 


is the winding number of the Hamiltonian, and c is a 
close loop with k varying from 0 to 27 r. The wind¬ 
ing number describes the total number of times that the 
Hamiltonian travels counterclockwise around the origin. 
Despite that the Zak phase only takes 0 or 1 with a 
modulus of 27 t, the winding number can be any integer, 
and indicates different topological phases. From Eq. ©• 
one can see that the definition of v is invalid whenever 
\h\ = yjh 2 + h 2 = 0 . As the energy spectrum of H{k ) 

is given by E(k) = h 0 ± ^Jh 2 + h 2 , it is obvious that the 
winding number is ill defined at the phase transition point 
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Fig. 1: (Color online) The sketch of the integral path of the 
winding number. 70 is the phase transition point, and the red 
dots stand for the degenerate points at which |/i(fc)| = 0. 

( 1 / 2 ) is the winding number of the system on the right (left) side 
of the dashed line with r] = r/o- The integral path is smoothly 
moved from one topological phase to another without passing 
the degenerate points. 


770, where the gap between two bands closes and one has 
h x (k 0 ,r]o) = h y (k 0l ri 0 ) = 0. 

To characterize the topological properties of the phase 
transition points, we define the winding number I'd on a 
circle around the gap closing point in the parameter space 
of k and 77. Following reference E3, we set k = A sin 6 + ko 
and 77 = A cos $+770, with 9 the varying angle, A the radius 
of the circle, and (k 0, 770) the gap closing point in the pa¬ 
rameter space of k and 77. Hence the Hamiltonian around 
the circle can be represented as h(9), and the winding 
number is defined as 



kixdkjy kjy(jtlTix 



( 6 ) 


with d a close loop from 9 = 0 to 9 = 2 -k. 
Setting h x /\h\ = cos a, one can easily have 



(7) 

( 8 ) 


which directly indicates the times that the Hamiltonian 
goes around the origin along the close loop c (or c'). Given 
a quantum state in the regime of 77 > 770, the correspond¬ 
ing winding number v\ is only defined on the momentum 
k space. As there are no singular points in the space of 
(k, 77) except of the gapless transition points schematically 
marked by dots in FigJT] one can analytically continue the 
integral path c to the two-dimensional parameter space.No 
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matter how we change the integral path, the winding num¬ 
ber is unchanged as long as the Hamiltonian does not cross 
the degenerate points, at which \h(k)\ =0. As schemati¬ 
cally displayed in Fig[l] one can smoothly move the inte¬ 
gral path pass the phase transition point. Such a contin¬ 
uous change of integral path does not change the value of 
the integral, but leaves close loops around each degenerate 
point and a close path only on the momentum space (as 
shown in Fig[T|(d)). Effectively, this process leads to 

y \vdj =Vl- V2, (9) 

i 

where v\ and v 2 denote winding numbers of different 
phases on two sides of the transition point, and de¬ 
notes the winding number of the i-th degenerate points. 
The above equation indicates clearly that the change of 
winding number v across the phase transition point equals 
to the summation of winding number v,i of each degener¬ 
ate points. 

Extended Kitaev model. — To give a concrete ex¬ 
ample and verify our scheme, first we consider the ex¬ 
tended version of the Kitaev model [20] by adding next- 
nearest neighboring hopping and pairing terms, descried 
by the Hamiltonian 

H = y tic|cj + i + t 2 c\c i+ 2 +h.c. 

i 

+AiCjCj + i + A 2 CiC i+ 2 + h.c. - 2fj,c\ci, (10) 

where cj (cj) is the creation (annihilation) operator of 
fermions at the i-th site. 

After the Fourier transformation, it takes the form of 

H = ^k H ( k )' l Pk, ( 11 ) 

where = (ck,c'_ k ) and 

"<*>=(& :?;)• <12) 

with h y = Aisinfc + A2sin2fc and h z = — ficosfc — 
t 2 cos 2 k + fj,. The eigenvalues are given by 

E(k)=±\h(k)\=±y/hl + hl (13) 

This model is related to the extended quantum Ising model 
with additional three-body interaction by Jordan-Wigner 
transformation mm- For simplicity, in the following dis¬ 
cussion we shall set Ai = t\ and A 2 = t 2 with t\ t2 taking 
real numbers [Tj. In general, according to the the ten¬ 
fold-way classification m , the model with complex A 12 
belongs to the D class, which shall be characterized by a 
Z 2 quantity. However, for the case with real Aj 2 con¬ 
sidered in the present work, due to the absence of a x in 
the Hamiltonian, this model belongs to the Z-type topo¬ 
logical system which can be characterized by the winding 
number. 


(b) 1 (c) (d) (d) 
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t 2 = -1.5 t 2 = -1 t 2 = -0.5 t 2 =0 



(f) (g) (h) (i) ' 
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t 2 = 0.5 t 2 =1 t 2 = 1.5 tj = 2 

Fig. 2: (Color online) The energy spectrum under OBC ver¬ 
sus t 2 and the winding of the Hamiltonian with different t 2 of 
the extended Kitaev model. The parameters of (a) are fj, = 1, 
ti = Ai = 0.5 and A 2 = t 2 . The times that the Hamilto¬ 
nian goes around the origin indicate the winding number of 
the corresponding topological phase. In (c), (f) and (h), the 
Hamiltonian goes through the origin, and the winding number 
is ill defined. 

Since the energy gap must close at phase transition 
points, we can determine the phase boundaries of the ex¬ 
tended Kitaev model by the gap closing condition |/i(/c)| = 
0, which yields (k = 0, ti+t 2 — fx = 0), ( k = 7 r, — ti+t 2 — 
/i = 0) or (k = ± arccos (— t\/2t 2 ), t 2 = —fx). Under open 
boundary condition (OBC), this model may hold different 
number (0, 1 or 2) of Majorana zero modes at each end. In 
FigEJa), we show the energy spectrum under OBC with 
/i = 1, t\ = Ai = 0.5 and A 2 =t 2 . Here we only display 
the upper branch of the spectrum with E > 0 due to the 
chiral symmetry of the spectrum. With increasing of t, 2 . 
we can see that the system goes through different phases 
labeled by the appearance and the disappearance of zero 
modes. Notice that there are two pairs of Majorana zero 
modes in the regime of t 2 < —1 or t 2 > 1.5, and only one 
pair in the regime of 0.5 <t 2 < 1.5. 

Next we calculate the topological invariant of this 
model, we rotate the Hamiltonian to the x-y plane with 
h y t h x , h z —> h y to coincide with the definition of the 
winding number in previous section. In Fig|2](b)-(i), we 
illustrate the winding of the Hamiltonian with different 
t 2 . When the momentum k varies from 0 to 27 t, the curve 
of the Hamiltonian may enclose the origin twice, once, or 
not enclose the origin, which correspond to v = 2, 1, 0, 
respectively. The number of times that the Hamiltonian 
goes around the origin indicates the number of Majorana 
zero modes at each end of an open chain. However, at the 
phase transition points (as shown in FigJ2|c), (f) and (h)), 
the curve crosses the origin, and the winding number is ill 
defined due to |h| = 0 at the origin. 

To characterize the topological properties of the phase 
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- v d — 1 - v d =—1 

- v d =-1, two degenerate points 



Fig. 3: (Color online) The phase diagram of the extended Ki- 
taev model with the winding number Vd on the phase bound¬ 
ary being marked. Here we have taken /z = 1, t\ = Ai, and 
£2 = A 2 . There are two degenerate points in the BZ for the 
system on the red dashed line, and the summation of Vd is 
—2. Inserts show the energy spectrum in momentum space for 
systems on the phase transition points labeled by stars. 



A = -2.5 




A = -1.5 


(d) 



A = -1 


(e) ; (f) ; (g) 




A = 0 A = 0.5 A = 1 


A = 1.5 


Fig. 4: (Color online) The winding of the Hamiltonian of the 
extended SSH model with different A and <5 = —0.5. The 
times that the Hamiltonian goes around the origin indicate the 
winding number of the corresponding topological phase. In 
(b), (d) and (g), the Hamiltonian goes through the origin, and 
the winding number is ill defined. 


When t' = 0, the model reduces to the SSH model [22] of 
the BDI class, which can be described by the Hamiltonian: 

H = 'Yytic\ i c B ,i + t 2 c\ i+l c B ,i] + h.c ., (17) 


transition points, we calculate the winding number v,i de¬ 
fined in the previous section by choosing £2 as the phase 
transition driving parameter, and demonstrate the phase 
diagram of this model with /z = 1 in Fig{3] We can see 
that the summation of Vd equals to the change of v across 
the phase transition point. The inserts (i), (ii) and (iii) in 
Fig{5]show the energy spectrum of phase transition points 
corresponding to Fig[2](c), (f) and (h), respectively. There 
is only one degenerate point with v& = 1 in (ii) and (iii), 
but two degenerate points emerge in (iii) with Vd = — 1 for 
each of them. While the change of v on two sides of the 
dashed-dotted lines is 1, the change of v across the dashed 
line is —2, which corresponds to the summation of v ( j of 
these two degenerate points. 


where c 


Extended SSH model. — Next we consider an ex¬ 
tended SSH model [H|, which takes the form of 


H = iplHikWk, 


(14) 


where i/>I = (ct a» c Lb) and 


t 

MB), 


is the creation operator of fermion on i-th 


H(k) = 


t\ + t 2 e lk + t'e 


t\ + t 2 e lk + t'e 


„-2 ik 


ik 


Alternatively, H{k ) can be expressed in the form H{k) = 
h x cr x + hyUy, with h x = t\ + t 2 cos k + t' cos 2k and h y = 
t 2 sin k + t' sin2fc. Diagonalizing the Hamiltonian, we get 
the eigenvalues 


E{k) = ±\h{k)\ = ±^hl + h? y . 


A (or B) sublattice. It is well known that this model has 
two topologically distinct phases. For convenience, we set 
ti = t(l—6), t 2 = f(l+<5) and take t = 1 as the energy unit. 
When S > 0, the system belongs to a topological nontrivial 
phase with the winding number v = 1, and supports two 
degenerate zero-mode edge states for the model with OBC. 
When 5 < 0, the system has v = 0 and the edge states 
disappear. 

When t' ^ 0, the system still belongs to the BDI class 
as the time reversal symmetry and the chiral symmetry 
are preserved, but the topological phases of the extended 
SSH model are enriched. Defining t' = At, the gap closing 
conditions are given by (k = 0, A = —2), (k = n, A = 25) 
or (fc = ± arccos [— (1 + <5)/2A], <5 + A = 1). By fixing 
5 = —0.5, we show the winding pattern of the Hamilto¬ 
nian with different parameters in Fig]4] It is clear that 
the winding pattern shows different geometrical structure 
in various parameter regimes. For the Fig£Ha) and (h), 
we can see that the curve goes around the origin twice, 
which corresponds to the winding number v = 2. Simi¬ 
larly, we have v = 1 corresponding to Fig]4][c) and v = 0 
corresponding to FigJ3](e) and (f). On the other hand, 
. (15) Fig0](b), (d) and (g) correspond to some phase transition 
points, at which the Hamiltonian goes through the origin 
and the winding number v is ill defined. 

Next we display the phase diagram of the extended 
SSH model in Fig[5] with topologically different phases 
characterized by different values of the winding number 
v. By choosing A as the driving parameter, we can cal¬ 
culate the winding number vy, for the phase transition 
points. As shown in Fig]5] we can see that v,i equals 


(16) 


P-4 





















Winding numbers of phase transition points for one-dimensional topological systems 


- v d —1 - v d =—1 

- V d =1, two degenerate points 
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Fig. 5: (Color online) The phase diagram of the extended SSH 
model with on the phase boundary being marked. There 
are two degenerate points in the BZ for the system on the red 
dashed line, and the summation of u d is 2. 

to the change of the winding number v across the phase 
transition point. The state at the red dashed line with 
5+A = 1 is a ID semimetal with band degenerate points at 
k = ± arccos [—(1 + 5)/ 2A], while the ones at other phase 
transition boundaries have only one degenerate point in 
the BZ. 

Summary. — In summary, we have studied the topo¬ 
logical properties of phase transition points of ID topo¬ 
logical superconductor and insulator systems. We charac¬ 
terized different topological phases with the winding num¬ 
ber, and defined the winding number for phase transition 
points in the parameter space of system’s momentum and 
a transition driving parameter. We studied several ID 
topological models with high winding number, and demon¬ 
strated that the topological phase transitions can be char¬ 
acterized by the introduced winding number around the 
transition point, the summation of which equals to the 
change of the winding numbers of different phases across 
the transition point. Our theory provides a way to classify 
different topological phase transitions by directly studying 
the properties of the phase transition point and can be ap¬ 
plied to study other ID Z-type topological systems. 
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